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Abstract: 

We have calculated the electroweak 0(a) corrections to the processes pp — ■> Z7 + X — > 
viva + X, pp -y Z7 + X -y U>y + X, and pp -> W7 + X -> uj^f + X at the LHC, with 
I = e, /x and ^ = z/ e , v T . The virtual corrections are evaluated in leading-pole approx- 
imation, whereas the real corrections are taken into account exactly. These corrections 
are implemented into a Monte Carlo generator which includes both phase-space slicing 
and subtraction to deal with soft and collinear singularities. We present numerical results 
for total cross sections as well as for experimentally interesting distributions. Applying 
typical LHC cuts, the electroweak corrections are of the order of —5% for the total cross 
sections and exceed —20% for observables dominated by high centre-of-mass energies of 
the partonic processes. 
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1 Introduction 



One of the primary goals of future colliders is to look for effects of theories beyond 
the Standard Model (SM). In general, there are two possibilities for new physics. Either 
the novel phenomena manifest themselves at energy scales probed directly at forthcoming 
experiments or at much larger scales. In the first case one could admire a spectacular 
scenario, characterized by the appearance of many new particles or resonances. In the 
latter case, the situation appears less favourable, but new-physics effects could still be 
detected in an indirect way. They can, for instance, influence physical observables at low 
energy by modifying the structure of the gauge-boson self-interactions. Such modifications 
can be parametrized in terms of anomalous couplings in the Yang-Mills vertices. While the 
gauge-boson couplings to the fermions have been measured at LEP2 and Tevatron to an 
accuracy of 0.1-1% [H], triple and quartic gauge-boson couplings have been determined to 
much lower accuracy [dl2]- Hence, the possibility of the existence of anomalous couplings 
in the electroweak gauge sector cannot be ruled out yet. 

Vector-boson pair-production processes turn out to be especially suited for testing the 
non-Abelian structure of the SM. In the last decade, their potential has been extensively 
exploited at LEP2 and Tevatron. For these processes, the effect of anomalous couplings 
is expected to rise strongly with increasing invariant mass of the gauge-boson pair, since 
non-standard terms lead in general to unitarity violation for longitudinal gauge bosons. 
At future colliders, it will thus be worthwhile to analyse the large di-boson invariant-mass 
region. Another interesting property of the vector-boson pair-production processes is the 
so-called radiation zero, which characterizes the SM partonic amplitudes of W7 and WZ 
production. Owing to gauge cancellations, for a given value of the scattering angle in the 
di-boson rest frame, the amplitude for W7 production vanishes exactly and the one for 
WZ approximately. For hadronic processes, the resulting effect is the appearance of a dip 
in the angular distributions (see e.g. Ref. and references therein). Any deviation from 
the SM gauge structure tends to modify this peculiar signature, thus providing a clean 
possibility to look for new physics. Present experiments are not sensitive enough to see 
the dip [HE], but the LHC offers the possibility to observe it for the first time. 

Finally, the production of vector-boson pairs constitutes an important background in 
the search for new particles. For example, 7Z production with the Z boson decaying into 
a pair of neutrinos gives rise to a photon-plus-missing-energy signature. The same signal 
is expected for processes where the photon is produced in association with one or more 
heavy particles which escape the detector, either because they are weakly interacting or 
because they dissociate into invisible decay products. Many extensions of the SM predict 
the existence of such processes. For example, the ADD model [ E] and related models [ 
|7J IB] , where gravity becomes strong at TeV energies, predict the possibility to observe the 
production of a photon along with a graviton. Higher-dimensional gravitons appear as 
massive spin-2 neutral particles, which cannot be directly observed in a detector. Hence, 
graviton radiation leads to a missing-energy signature where a photon is produced with no 
observable particle balancing its transverse momentum. Independently of the considered 
new-physics model, the production of such heavy invisible particles can be inferred from 
the missing transverse energy distribution, which would exhibit an excess of events at 
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high values as compared to the SM. Existing experimental searches for such an effect are, 
for instance, presented in Refs. ITU]. 

Owing to their limited energy and luminosity, LEP2 and Tevatron only provide weak 
constraints on anomalous couplings and production rates for new particles [IT1 ITT1IT21 IT*"] . 
At the LHC, the experimental collaborations will collect hundred thousands of events 
coming from vector-boson pair production [ [*["""| • To match the precision of the LHC 
experiments, the cross sections for vector-boson pair-production processes have to be 
calculated beyond leading order. The next-to-leading-order (NLO) QCD corrections for 
on-shell W and Z bosons were calculated in Refs. [ El EH) and extended to include leptonic 
decays in the narrow-width approximation and anomalous couplings in Refs. [ I17| IT"""! [T§]. 
After the full NLO amplitudes including leptonic decays had become available [ l2T"| . 
Monte Carlo programs incorporating these amplitudes have been presented for W7 and 
Z7 production in Ref. [EU and for WW, WZ, and ZZ production in Refs. [I221I23]- As a 
general feature, the QCD corrections modify the leading-order di-boson cross sections at 
the LHC by a positive amount of several tens of per cent, and have thus a considerable 
impact on the measurement of gauge-boson couplings. At this point, the question arises 
if the electroweak corrections to vector-boson pair-production processes at the LHC lead 
to similarly sizable effects. In the above-mentioned calculations typically only universal 
electroweak corrections were taken into account such as the running of a or corrections 
associated with the p-parameter. This approach was based on the idea that nonuniversal 
electroweak corrections only contribute significantly in the high-energy domain, where the 
statistics at the LHC will be too low to extract any physics from the experimental data. 

However, this is in general not the case in the high-energy region, which is of consider- 
able experimental interest since as discussed above the effects of the anomalous couplings 
and the effects of the production of so far unknown heavy particles are most pronounced 
there. In fact, at high energies the electroweak corrections are enhanced by so-called 
electroweak Sudakov logarithms, i.e. double logarithms of the process energy over the 
vector-boson mass ["2H 123 I2U1 12*71 I28| . They can reach several tens of per cent and have 
to be taken into account to make sure that an experimentally observed deviation from 
the QCD-corrected SM predictions due to electroweak effects is not misinterpreted as a 
signal for physics beyond the SM. The importance of the electroweak corrections to gauge- 
boson pair-production processes at the LHC has been confirmed by several calculations 
[ [2111 EII1 H*"2 > which all point to a negative contribution of some ten per cent, which well 
exceeds the statistical errors. Whereas the logarithmic electroweak corrections to WW, 
WZ, ZZ production in the high-energy limit have been calculated including the corrections 
to the decay process as well as the complete real corrections [EJ, f° r W7 production only 
the logarithmic electroweak corrections to the production subprocess are available [ l2*"J]. 
The electroweak corrections to Z7 production have so far only been determined for an 
on-shell Z boson, i.e. without taking into account decay effects [|3U]. Furthermore, the real 
corrections for the latter two processes have been included only in the soft and collinear 
limits or not at all. The goal of our work is therefore the calculation of the virtual and 
real electroweak corrections to the cross sections and distributions for the purely leptonic 
channels pp — > + X — > li^ + X, where V = W, Z. Final states containing quarks are 
not discussed in this work. 
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This paper is organized as follows: The strategy of our calculation is described in 
Section |21 In Section El we give the setup for the numerical evaluation. Numerical results 
are presented in Section and Section contains our summary. Finally, some explicit 
analytical results are listed in the appendices. 

2 Strategy of the calculation 

We consider the production of a photon along with a massive gauge boson in proton- 
proton collisions, where the gauge boson decays leptonically. The generic process reads 

p + p^y + 7 + X^/ 3 + r 4 + 7 + X, (2.1) 

where p denotes the incoming protons, V indicates a W or Z boson, 7 the outgoing photon, 
l 3 the outgoing lepton, l 4 the outgoing anti-lepton, and X the remnants of the protons. 

In the parton model, the corresponding cross section is obtained as a convolution of 
the distributions & qi \ p and $g 2 |p °f ^ ne partons q 1 and q 2 in the incoming protons with 
the partonic cross sections <T qiq2 , averaged over spins and colours of the partons, 

cr pp( 5 ') = / dxi / dx 2$ qi |p(xi)$5 2 |p(£2)cw(s)- (2.2) 
11,12 Jo Jo 

The sum Yl qi q2 runs over an partonic initial states which are allowed by charge conser- 
vation. In practice we include the quarks and antiquarks g« = u, u, d, d, s, s, c, c. The 
integration variables Xi,i = 1,2 correspond to the partonic energy fractions. Moreover, 
the quantities S and s = X1X2S are the squared CM energies of the hadronic process and 
the partonic subprocesses, respectively. 

For the partonic processes we adopt the generic notation 

<?i0i,ci) + 92(^2,0-2) -> V(p 3 +p 4 ) +7(P5,A 5 ) -> l 3 (p 3 , <r 3 ) + [4(^4,^4) +7(P5,A 5 ), (2.3) 

where the arguments label the momenta and helicities <7j = ±1/2, i = 1,...,4 and 
A5 = ±1 of the corresponding particles. In the following, the masses of the external 
fermions are denoted by rrii, i — 1, . . . , 4. Since we treat all external fermions as massless 
these appear only as regulators of mass singularities. The charges of the fermions are 
denoted by Qi, those of the anti-fermions by —Qi- 

The lowest-order partonic cross sections are calculated using the complete matrix 
elements for the process 

?i(Pi, o-i ) + Q2(?2, cr 2 ) -> h(p 3 , cr 3 ) + r 4 (p 4 , cr 4 ) + 7(p 5 , A 5 ). (2.4) 

This means that we include the full set of Feynman diagrams, in this way accounting for 
all irreducible background coming from non-resonant contributions. Analytic results are 
given in App. IA.1I 

The electroweak radiative corrections to ()2.4j) consist of virtual corrections, resulting 
from loop diagrams, as well as of real corrections, originating from the process 

?i(Pi, + q 2 (p2, cr 2 ) -> h(p3, cr 3 ) + l 4 (p 4 , cr 4 ) + 7(^5, A 5 ) + 7(p 6 , A 6 ) (2.5) 
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with an additional photon with momentum pq and helicity A6 = ±1. Virtual and real 
corrections have to be combined properly in order to ensure the cancellations of soft and 
collinear singularities (cf. Section I2T2"|I . The helicity matrix elements for the process (J2.5I) 
are listed in App. IA.31 For the calculation of the radiative corrections we follow the 
approach used for the process e + e~ — > W + W~ — > 4/ in Ref. [E2] : The virtual corrections 
are calculated in the Leading-Pole Approximation (LPA), i.e. we take only those terms into 
account that are enhanced by a resonant massive gauge-boson propagator (cf. Section liOj) . 
whereas the real corrections are calculated from the full matrix elements for the processes 



2.1 Absorption of initial-state mass singularities in parton distribution func- 
tions 

At next-to-leading order in electroweak perturbation theory, the partonic cross sections 
generally contain universal, initial-state collinear singularities. The latter can be absorbed 
in the Parton Distribution Functions (PDFs). To this end, one has to fix a factorization 
scheme, and we choose the MS scheme, which defines the connection between the lowest- 
order PDFs in (|2.2|) and the experimentally determined PDFs as [ 031 1 



-srMS 



1 dz —MS (Xi 



a 



(2.6) 



Here, we used the (+)-prescription 



dx [f(x)],g(x) 



dxf(x)[g(x)-g(l)]-g(l) / dxf{x), 



(2.7) 



and the quantity Q denotes the electromagnetic factorization scale, which should be set 
equal to the typical energy scale of the considered process. Inserting the MS definition 
(12.61) into the hadronic cross section (|2.2j) allows to write the latter as a convolution of 
the experimentally defined PDFs with the subtracted partonic cross section 
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a pp (S,Q 2 ) = J2 I d*! / dx 2 $^ p (x 1? g 2 )$^ p (x 2 ,Q 2 X|(S,g 2 ), (2.8 



where 
and 



a 



a 9l9- 2 ( S ' Q ) = ° 
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qi 52, sing 



f£) - rr MS 

9192 V" 5 J u q 1 q 2 , sing 



[s,Q 2 ) + 0(a 2 



a 
4ns 



['dz [ d$ 4 (z) M® b (H*)) 

i=l,2 J ° J 



1 

X - 

z 



i±^(ln(g)-21n(l-,)-l 



(2.9) 



(2.10) 
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Here M qi q 2 denotes the lowest-order matrix element, and <&i(z) represents the partonic 
phase-space where the incoming particle i has momentum zpi. After the subtraction, the 
resulting cross section is free of collinear and soft singularities. 

A fully consistent inclusion of the electroweak corrections into the calculation of the 
hadronic cross sections requires the use of PDFs in which both photonic initial states and 
QED 0(a) corrections are taken into account. A first global analysis of parton distribu- 
tions incorporating QED contributions has been performed recently in Ref. [EH]- It turns 
out that the photon PDFs are quite small [EH], and the QED corrections change the PDFs 
of the quarks by less than 1% [EH1 EZj - These effects are below the typical uncertainty 
of hadronic processes and below the envisaged accuracy of our analysis. Accordingly, we 
restrict our calculation to initial states containing quarks and use only QCD-corrected 
PDFs in (12. 8|) . We nevertheless investigated the dependence of our results on the electro- 
magnetic factorization scale Q entering in ()2.10|) . keeping the QCD factorization scale Q 

that enters $^|p(xj, Q 2 ) in (J2.8|) fixed. For 7Z production, a variation of Q by a factor 10 
resulted in a change of the total cross section by less than 0.3%, of the 7Z invariant-mass 
distribution at large invariant mass by less than 0.5%, and of the transverse momentum 
distribution of the photon at large transverse momenta by less than 1%. 

2.2 Treatment of soft and collinear photon emission 

In first-order electroweak perturbation theory, the subtracted partonic cross section 
receives contributions from the lowest order, the virtual corrections, the real corrections, 
as well as of the term (j2.10|) arising from the factorization of the initial-state collinear 
singularities: 

<|(*> = + + &) - <Ling(*, Q 2 )- (2-11) 

In general, virtual and real corrections contain soft as well as collinear singularities, which 
we regularize by introducing an infinitesimal photon mass m 7 as well as small mass param- 
eters rrii for the external fermions. These logarithmic singularities cancel when adding 
real and virtual corrections, except for the collinear logarithms arising from the initial 
state, which are absorbed in the PDFs as described in Section |2~T1 Thus, all soft and 
collinear singularities cancel in (j2.11|) . leaving the subtracted cross section free of any 
mass regulators. 

However, this cancellation requires a careful treatment of soft and collinear singular- 
ities. In our work, we adopt the approach to divide the virtual and real corrections into 
finite and singular parts, and to cancel the dependence on the mass regulators by adding 
the singular parts analytically. To extract the singular parts from the real corrections, 
we employ two independent methods, the phase-space-slicing technique and the dipole 
subtraction method. In the following we give explicit expressions for our treatment of 
infrared, i.e. soft and collinear, singularities. 



2.2.1 Finite virtual corrections 

First, we define the finite virtual corrections by 



virt ( ~\ 

°<3i<?2, finitely 



2l ' ' 
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-2 Re 



(2.12) 



where A4^g 2 and denote the matrix elements of the Born process and the corre- 

sponding electroweak virtual corrections, respectively, whereas <3> is the phase-space for 
the Born process ([2.4)1 . The subtracted singular parts are given by 



2 Re 



a 
27 



(2.13) 

where n is the number of external particles of the lowest-order process, s^- = (Pi+Pj) 2 , and 
Tj = 1 for incoming particles and outgoing anti-particles, whereas r« = — 1 for outgoing 
particles and incoming anti-particles. The singularities are contained in 



C(sij,m 2 ) 



ml 
In I 



1 + lnf^ 



2 \ Sy 



ilnf^ 

2 y Sjj 



;2.i4) 



while the finite parts are added for later convenience and distinguish between incoming, 
z, j = 1,2, and outgoing particles, z, j > 2, 



7T 



7T 



7T 



3^2' 

2 3 
h - 

2^2' 

i + 2, 



i = l,2, J>2, 
i>2, J = 1,2, 
i = 1,2, J = 1,2, i^j. 



(2.15) 



£.£.,2 Phase-space slicing 

In short, the idea of the phase-space-slicing technique is to define the finite real correc- 
tions by restricting the phase-space integration to the region where the squared amplitude 
of the bremsstrahlung process is finite. To this end, a technical cut 5 S on the energy E 
of the bremsstrahlung particle is introduced, as well as a technical cut 5 C on the angles 
di of the bremsstrahlung particles with respect to all potential emitters. The finite real 
corrections are then given by 

<s:l.« = ± I , ^ d *- i*s<*->f. (2.i6) 

cos i9 j < 1 — S c 

where $ rea i and .M real denote the phase- space and the matrix element of the brems- 
strahlung process (|2.5J) . 

In the remaining singular parts of the real corrections, the momentum of the brems- 
strahlung particle is integrated out, up to a remaining convolution over the momentum 
fraction z of the initial state particle after bremsstrahlung emission. Here the soft and 
collinear singularities are again regularized by an infinitesimal photon mass m 7 and by 
small mass parameters for the external fermions. The resulting expressions 
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can be divided into a part arising from initial-state radiation and a part arising from final- 
state radiation. Adding them to the term (|2.1(Jj) originating from the MS definition of the 
PDFs and to the singular parts (|2.12j) of the virtual corrections, all infrared singularities 
cancel and one obtains a finite contribution 



virt+real,sli / ^2 
°9ig2,sing l S ' V 



virt / £\ , real,sli / ~\ _ -.MS , 
°gi92,singV 6 / ~r °9i92,singV 6 / a q 1 q 2 ,smg' 



s,Q 2 



2s 



d$ r 



^v+r.init 



9ig 2 ,singV - 
2 



'gi 92, sing \ 



i=l ^ ^ 



gi 92, sing 



, (2-17) 



where 



9ig2,singV^U / 



2 n 



i=l 3=1 

35"^ 



X 



7T 2 _ / AE,E 

— - 2 + dj ; + Li 2 1 - 



+ ln 



i 



2 . ln«l+|ln« 



(2.18) 



A^v+r, final 
J l c/i92,sing 
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27 



i=3 3=1 

35"^ 



.M (0) 



'/.I </■-> 



X 



7T 



- + Qj + Li 2 



+ ln ' 



4^- 
1 



3 i 
- + In , 

2 V 4£ 2 



-In 2 



4£; 



(2.19) 



'9i 92, sing 



x 



3=1 



2z 



^( ln (_i| j + 21n(1 _ z) + 1/ 



(2.20) 



Li 2 denotes the dilogarithm 



dt 

Li 2 (z) = - I — ln(l-t), | arc(l - z)\ < ir, 



(2.21) 



and 3>i(#) represents the phase-space where the incoming particle i has momentum zpi 
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2.2.3 Dipole subtraction method 



In the dipole subtraction method, the finite real corrections are constructed by sub- 
tracting an auxiliary function from the squared bremsstrahlung amplitude before inte- 
grating over phase-space, 



real, sub /~\ 
^gi g 2 , finite I S J 



2i 1 d$rea 



(2.22) 



The subtracted terms are added again after partial analytic integration over the brems- 
strahlung momentum. The subtraction function has to be chosen such that it cancels 
all soft and collinear singularities of the original integrand, so that the difference ()2.22|) 
can be integrated numerically without regulators for these singularities. Moreover, it has 
to be simple enough so that it can be integrated analytically over the singular regions 
of phase space. We use the process-independent dipole subtraction formalism, which was 
introduced for massless QCD in Refs. [ EHl Ell E] an d generalized to massive fermions in 
Ref. [EE]- We follow the approach of Ref. [EH] i n the limit of small fermion masses. 

In the dipole subtraction formalism the subtraction function is constructed from con- 
tributions that are labelled by ordered pairs ij of charged fermions, so-called dipoles. 
The fermions i and j are called emitter and spectator, respectively. In the formulation of 
Ref. [EHl the subtraction function reads 



l^sub 



'11 >J2 



real J 



Pj i Pn b 



= 5 !J = 1 



M {0) - CC fc 



(2.23) 



The sum over n& = 5,6 accounts for the fact that the photon appears twice in the final 
state of the bremsstrahlung process (|2.5|) . such that a separate subtraction term has to be 
introduced for each final-state photon. Note that one of the two final-state photons 
has always to be visible in the detector and thus does not give rise to singularities. 
In what follows, we suppress the dependence of the phase-space on n& and denote the 
bremsstrahlung momentum by k = p n , b . The phase-spaces i n Q2.23j) are given by 
embedding prescriptions of the phase-space $ real of the bremsstrahlung process in the 
phase-space $ of the Born process. Both, the embedding prescription • and the 
subtraction functions g™ b depend on the emitter % and spectator j being incoming or 
outgoing. One therefore encounters four kinematically different cases for the subtraction 
functions gf" b 



(pikXl-yij) 
1 



i 
i 

, (Pik)vij 



l-Zij(l-yij) 
2 



2-x 



2 



X 



Ji 



1-Vi 



i,j > 2, 

i>2, j = 1,2, 
i = l,2, j>2, 
i = l,2, j = l,2. 



(2.24) 
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where 



J. ; 



PiPj +Pjk — pjk 
PiPj + Pjk 
PiPj 



Pik 



PiPj + Pik + Pjk 
PiPj - Pik- Pjk 



PiPj +Pjk' 



PiPj 



(2.25) 



Denoting the number of external particles by n, one has for the embedding prescription 
if both the emitter and the spectator are part of the final state 



l - y y J 



-1-1 

Pi 



1 



-Pi 



~u 

Pi 



1 7 ' 

i,3>2, i^J, 



Pi 



l = l,...,n, l^ij, (2.26) 



whereas for an initial-state emitter and a final-state spectator one gets 



Pi = X jiPi, Pj = Pj + k>1 ~ (1 - Xji)j%, 



Pi 



Pi 



1,2, j>2, 1 = 1,.. .,n, (2.27) 



The embedding prescription for a final-state emitter and an initial-state spectator can be 
obtained from (|2.27|) by exchanging % <-> j everywhere. Finally, in case of an initial-state 
emitter and an initial-state spectator the embedding prescription is given by 



Pi = V ijPi, Pj = Pj, 

with the boost matrix 

= 



Pi 



9\ 



= A"X, i,j = 1,2, i^j, I = 3,...,n (2.28) 

(Pij + Pij)^{Pij + Pij)u ^PijPij^ 



P.2 i p.. p. 



p.2 



(2.29) 



Here, P^ represents the total initial-state momentum after the projection, — pf + pj, 
and 

P lJ =Pi+P2-k = Y,Pk- (2-30) 

k=3 

The subtracted contribution can be integrated over the (singular) photonic degrees of 
freedom up to a remaining convolution over a variable z. In this integration the regulators 
m 7 and m ; must be retained, and soft and collinear singularities appear as logarithms 
in these mass regulators. The resulting expressions [ E21 EE] can again be divided into 
a part arising from initial-state radiation and a part arising from final-state radiation. 
Adding them to the term (J2.10|) originating from the MS definition of the PDFs and to 
the singular parts (J2.12)) of the virtual corrections, the infrared singularities cancel and 
one finds 
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(2-31) 



9 



where 



'gi(j2,sing 



2 QiTiQjTj 

j'=l 
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a<£ 2 (<^)) 

+ 21n(l ) +1-2 



1 - z 



In ' l%l 



Q 2 z 
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, (2.32) 



v+r , final, z, 
qi qs, sing ' 



2^ 



1-z 



2 In 



(2.33) 



and <5j>2 = 1 for j > 2 and 5j>2 = for j < 2. Note that the finite part ()2.15|) of the 
singular virtual corrections ()2.13|) has been chosen to exactly cancel the end-point parts 
resulting from the subtraction function (|2. 23)1 such that only z-dependent contributions 
remain in (jZHISD and 



Master formula 

Altogether, in our approach the subtracted partonic cross section, which is free of soft 
and collinear singularities, can be written as a sum of four different parts 

= "fiM + <«(*) + <lfin ite (^) + *££K(*' (2-34) 

where the finite virtual corrections in the second term of the right-hand side are defined 
by (J2.12|) . The last two terms are given by ()2.16|) and (|2.17|) . respectively, if using the 
phase-space-slicing technique and by (|2.22|) and (j2.31|) . respectively, in case of using the 
subtraction method. 



2.3 Treatment of finite- width effects 

The processes (|2.3|) involve the production and decay of an unstable particle V = W, Z. 
The corresponding propagator leads to a pole in the amplitude, which has to be regularized 
by incorporating the finite width of the massive vector bosons. Several methods to include 
finite-width effects in a perturbative calculation have been discussed in the literature [ 
g2lE2J EH US Eni 113 IIHI SHI HOI EH H21 EH] The simplest approach is given by the fixed- 
width scheme, which corresponds to replacing the resonant propagator by one containing a 
constant finite width. Since the finite width results from resummation of an incomplete set 
of higher-order contributions, its introduction potentially violates gauge invariance. For 
Born amplitudes, the effects of this violation have been shown to be numerically small in 
the fixed-width scheme (see e.g. Refs. [HHlllZlEllEn]) by comparing the results to those 
obtained with other manifestly gauge-invariant prescriptions, as e.g. the complex-mass 
scheme [HZ]. 
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In our calculation, following the approach of Ref. [|o2], we adopt the fixed-width 
scheme to evaluate the Born cross section and the real corrections. For the virtual cor- 
rections, on the other hand, we use the LPA [ |321 ESI EH] • In this approximation only the 
leading term in an expansion about the resonance pole, i.e. the residue divided by the 
resonant propagator, is kept. Since the residue, which is related to physical amplitudes 
for on-shell production and decay, and the resonant propagator are gauge-invariant, the 
LPA yields a gauge- invariant result. As all non-resonant contributions are neglected, the 
number of contributing diagrams is considerably reduced, and only loop integrals with up 
to four propagators appear. The error induced by using the LPA for the virtual correc- 
tions can be estimated to be of the order aTy/My, where My and Ty are the mass and 
the width of the vector boson V, respectively, and thus of the order of a few per mille. 
This estimate was confirmed for single W production at hadron colliders in Ref. [ 02| and 
for WW production in electron positron annihilation in Refs. [ EEl EZ] • Thus, the LPA 
should be more than sufficient for the considered process, as long as the resonant diagrams 
dominate. 

Following the approach presented in Ref. [E2]; the virtual corrections evaluated in 
LPA can be further split into factorizable and non-factorizable parts, such that the finite 
virtual corrections can be written as 

<^ZM = ^/d$ 2Re [<r>o,$o sh ) 

x (&<££ A (*o, $o sh ) + SM^i(^ , $r ) - SM^ S (%, $o sh ))*] • (2-35) 

In f[2~35|) . A^ 2 LPA ($ , $g sh ) denotes the Born matrix element in LPA 

A/f(°)' LPA f5) ^oshN R{0) ( $ o sh ) /o oe\ 

M ™ ($ °' $0 >~ fa+ptf-Mt+iMyVv' (2 ' 36) 

which receives contributions from all Born diagrams containing the resonant propagator. 
The residue R of the amplitude at the pole is evaluated using a phase-space $Q Sh = 
{p° sh , i — 1, ... ,5} projected on the mass-shell of the decaying particle, i.e. (p3 Sh +p4 Sh ) 2 = 
M v . We specify our choice for the on-shell projection in App. [Bj The factorizable virtual 
corrections in ()2.35|) get contributions from all diagrams where the resonant propagator 
appears outside the loop, i.e. from diagrams that factorize into a production part, a decay 
part, and the resonant propagator. In LPA, they can be written as 



pvirt/^osh m 

*M Mc ($ , $ ) - {p3+PA y-M v + iMyr v > (2 ' 37) 



where all soft singularities are regularized by the infinitesimal photon mass parameter m 7 . 

The non-factorizable virtual corrections receive contributions from all diagrams where 
a massless virtual particle, in our case a photon, connects the production and decay 
subprocesses or one of these subprocesses to the resonance. The infrared singularities 
originating from photon emission off the unstable boson V are here regularized by keep- 
ing the resonant momentum off-shell wherever the on-shell limit leads to a singularity. 
The contribution of the diagrams where the photon couples to the resonance are already 
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partly contained in the factorizable virtual corrections as defined above, such that this 
contribution has to be subtracted from the non-factorizable corrections to avoid double 
counting [E2]- As an advantage, the non-factorizable virtual corrections can be evaluated 
in the Extended Soft Photon Approximation (ESPA), which corresponds to setting the 
loop momentum to zero wherever this does not lead to a resonance. The non-factorizable 
corrections 

^• nfacl °' ' (p 3 + Pa) 2 ~M v + iM v T v 



2 



M^ PA ($ , n sh )&K\ (pa + P.) 2 ~ My), (2.38) 



are then proportional to the Born matrix element in LPA. We give the analytical expres- 
sion for the correction factor 5™^ c in App. IA.2I 

Note that in order not to spoil the cancellation of soft and collinear singularities 
between the virtual corrections in LPA and the real corrections evaluated exactly, we 
consistently subtract in ()2.35|) the singular virtual corrections evaluated in LPA 



2 Re 



2 



f- ^nrjQiQi ; MS 2 LPA (^o,^ sh ) (jOis^m^ + Qj), (2.39) 



2tt 



osh 
ij 



but we use the exact singular virtual corrections in (|2.17|) and ()2.31jl . In ()2.39|) . s 
(p° sh + p° sh ) 2 , and £ and Cy are given by (J2.14j) and (|2.15|l . respectively. 

2.4 Implementation 

We have implemented our strategy in a Mathematica package called Pole, which works 
as an extension of the computer-algebra packages FeynArts3 [EH] and FormCalc3. 1 [ 155] . 
To be explicit, Pole extends these packages by the following features: 

• amplitude generation for initial-state hadrons, 

• amplitude generation for a given intermediate state to select only resonant factor- 
izable and resonant non-factorizable diagrams, 

• evaluation of the non-factorizable part of the amplitude in ESPA 

• reduction of fermion chains using the Weyl-van der Waerden formalism, 

• numerical evaluation of hadronic cross sections, 

• numerical evaluation of the virtual corrections in LPA, 

• phase-space integration for an arbitrary number of external particles using the 
generic phase-space generator of Lusif er [IfiO]. 
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• numerical evaluation of the real corrections using phase-space slicing or the dipole 
subtraction method. 

The main changes for the amplitude generation are the possibility to specify hadronic 
initial states as well as the resonant particles in the intermediate state. The latter infor- 
mation is used to only generate the resonant factorizable and the resonant non-factorizable 
diagrams. In the FormCalc part, we have implemented the Weyl-van der Waerden for- 
malism, which reduces the fermionic spinor structures down to Weyl-spinor products 
and Kronecker deltas in the helicities. Note that our implementation of the Weyl-van 
der Waerden method works independently of the formalism contained in FormCalc4 [ 
I5T] . where the fermionic structures are numerically evaluated as two-dimensional matrix- 
vector products. As an advantage of our approach with regard to FormCalc4, the matrix 
element is split into helicity amplitudes, and only non-vanishing components are calcu- 
lated. 

As a consequence, the kinematical dependence of the amplitudes as calculated by Pole 
is completely different from the amplitude given back by FormCalc3. 1, so that we had to 
rewrite the Fortran part of FormCalc3 . 1. The input required by the new code are model 
parameters, kinematical cuts, PDF set, as well as the aforementioned on-shell projection of 
the outgoing momenta. Using this input, the program automatically calculates hadronic 
cross sections and kinematical distributions including real QED corrections as well as 
virtual electroweak corrections evaluated in LPA. The phase-space integration is done 
using the generic Monte Carlo generator of Lusifer [ IHOJ , which is able to handle the 
kinematics of a process with an arbitrary number of final-state particles. To flatten the 
propagator peaks, Lusifer employs the multi-channel importance sampling [ 1521 153*1 I54"t 
153] with adaptive weight optimization [155]. Since the generator was originally designed 
to calculate Born cross sections only, we had to extend it to incorporate an appropriate 
treatment of soft and collinear singularities. We implemented both, the dipole subtraction 
method and the phase-space-slicing technique as presented in Section "2~!""l 

Pole is designed to handle processes involving an arbitrary final state and an arbitrary 
number of resonances. The applicability is limited to processes where the decay products 
are stable, i.e. it cannot handle cascade decays where unstable particles decay into unstable 
particles. Furthermore, the treatment of the non-factorizable corrections as well as the 
implementation of the phase-space slicing and subtraction methods are limited to massless 
external particles. 

3 Input and checks 

3.1 Parameter input and process definition 

We consider three classes of processes: 

p + p— >ui + Pi + j + X, l = e,/i,r, (3.1) 
p + p^Z+T+7 + A, / = e,/i, (3.2) 
p + p^^+T+7 + A, / = e, pL. (3.3) 

The first two classes allow to analyse Z7 production, while the third one contains W7 
as intermediate state. In the following sections, we present results for the LHC at CM 
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energy v^S* = 14TeV and an integrated luminosity L = 100 fb -1 per experiment. We 
neglect all fermion masses except for the mass of the top quark. We though keep the 
fermion masses in arguments of loop integrals as a regularization parameter for possible 
collinear divergences, where we use the values [ |H3 EE] 

m u = 0.066 GeV, m c = 1.6 GeV, m t = 178 GeV, 

m d = 0.066 GeV, m s = 0.15 GeV, m b = 4.9 GeV, 

m e = 5.109989 x 10" 4 GeV, m M = 0.105658369 GeV, m T = 1.77699 GeV. (3.4) 

For the vector-boson masses and decay widths we take [ IHZj 

M w = 80.425 GeV, T w = 2.124 GeV, 

M z = 91.1876 GeV, T z = 2.4952 GeV, (3.5) 

and the Higgs-boson mass is fixed to 

M H = 115GeV. (3.6) 

We use the G M scheme, i.e. we define the fine structure constant in terms of the Fermi 
constant 

« = « G , = ^G M M^4, (3.7) 

where the quantity = 1 — M^/M| denotes the sine of the weak mixing angle squared, 
and we use the value = 1.16637 x 10~ 5 GeV -2 for the Fermi constant. The definition 
(13. 7|) effectively resums contributions associated with the evolution of a to the W-boson 
mass and incorporates leading universal m t -dependent two-loop corrections. For the cou- 
pling of one real final-state photon, we use instead the value a = a(0) = 1/137.035999, 
i.e. we rescale the cross section with a(0)/aa „. Thus we have for the lowest-order cross 
section and for the corresponding corrections er^: 

a {0) oc a(0)c4 M , a (1) oc a(0)a 3 G ^. (3.8) 

The terms absorbed in «c M have to be subtracted from the finite virtual corrections ()2.35|) 
calculated in the a(0) scheme, once for each vertex parametrized by ac^, to avoid double 
counting. Since in our case we have two such couplings in the lowest-order cross section, 
the finite virtual corrections are evaluated according to 

<SSe, aGM (^) = - 2a^ P \s)Ar^. (3.9) 

The one-loop contribution Ar^ to Ar can for instance be found in Ref. [169]. 

We neglect all loop corrections to the quark mixing, and simply multiply the cross 
sections of the partonic subprocesses with the squares of the quark- mixing matrix elements 



Kd = 0.974, Ks = \ 1 - K 2 H , V nh = 0, 



Kd = -yi-^ d , Kb = 0.974, Kb = 0, 

Kd = 0, V tc = 0, V th = l (3.10) 
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for W7-production processes 

Following our discussion at the end of Section 12.11 we neglect the QED corrections to 
the PDFs. We therefore choose the only QCD-corrected CTEQ6M set [ I70j to calculate 
the hadronic cross sections and we neglect all contributions of the initial states containing 
bottom quarks. Denoting the momenta of the incoming protons by P\ >2 and the final-state 
momenta collectively by pfi n , the hadronic cross section for the Z7-production processes 
()3.1j) and ()3.2|) can be written as 



v PP (S,Q 2 ) = J d Xl J dx 2 2x ^ s J d$(p fi „) 

E (\Mgg(x 1 P 1 ,X 2 P2,P&n)\ 2 ^p(Xl, Q 2 )$^p(x 2 , Q 2 



X 

g=u,d,c,s 



+ \M qq (x 2 P 2 ,x 1 P 1 , Pta )\ 2 ^2( x ^Q 2 )^p( x ^Q 2 ))^ ( 3 - n ) 



and we use Ai C c = -M uQ and Ai sS = M-dA in the zero mass limit for the external fermions. 
For the W7-production process (j3.3|) . the effects of the quark mixing are taken into 
account by setting the quark-mixing matrix to the unit matrix for the calculation of 
the amplitudes, and by convoluting the squared amplitudes according to 



a pp {S,Q 2 ) = J dxi J dx 2 2x ^ s J d$( Pfin ) 



X 

g=u,c g'=d,s 



\M w (x 2 P 2 ,x 1 P 1 ,p &n )\^ q{ ;(x 2 ,Q^ ql]p (x u Q')) , (3.12) 



and we use -M ud " = Mus = M- C d = M.cs for unit quark-mixing matrix. 

Finally, for the typical energy scale Q used to evaluate the PDFs we choose [ HI] 

Q 2 = l -{M 2 v+P 2 Ty + P 2 ^) , \/ = Z,W, (3.13) 
where v denotes the transverse momentum of the resonance 



Pt,v 



Pt 



miss 



(PT, X +P l T, X ) 2 + (PT,y+PT,y) 2 ^1 PP iTy, (3.14) 

(p l T , x +P^) 2 + (PT,y+PT 



miss \ 2 



for 


PP " 




for 


PP " 


- ih, 


for 


PP " 


-> i>iij. 



The missing transverse momentum p™ lss is here defined as the sum of the transverse 
momenta of all particles visible in the detector, whereas p mi ^ s and p mi ^ s denote missing 
momentum in x and y direction, respectively. We use the same scale Q in the calculation 
of the electromagnetic MS subtraction term (|2.1()j) . 

We have implemented a general set of cuts, proper for LHC analyses. In order to 
specify this set, we introduce the azimuthal-pseudorapidity separation between particle i 
and j 



&R ij = J(T H -Vj) 2 + (<Pi-<Pi) 2 , (3-15) 
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where <fi is the azimuthal angle parametrizing the spatial part of the external momentum 
Pi, and the corresponding pseudo-rapidity 



rji — — In I tan ) (3.16) 



is a measure of the production angle 6i of the considered particle with respect to the 
beam. 

In order to define observables that are free of soft and collinear singularities, we use 
the following recombination scheme: 

• Photons close to the beam, i.e. with a rapidity |?7 7 | > ?7 rcc = 2.5, or with too small 
energy, i.e. E 1 < E rcc = 2 GeV, are treated as invisible, i.e. they only contribute to 
the missing momentum. 

• Photons with |?7 7 | < r] vcc = 2.5 and E 1 > E rcc = 2 GeV which are close to a final 
state charged lepton or another photon, i.e. which fulfil Ai? 7i < R rec = 0.1 for 
at least one i = e, /i, 7, are recombined. This means that the momentum of the 
photon and the momentum pi of the final-state charged lepton or photon for which 
AR^i has the lowest value are added and considered as an effective lepton or photon 
momentum. 

• Photons with |rf | < rf cc = 2.5, E 1 > E TCC = 2 GeV, and AR^ > R TCC = 0.1 for all 
i = e, /i, 7 are treated as visible. 

For the acceptance cuts, we consider an event to contribute to the cross section if the 
momenta for the leptons after possible photon recombination and the momentum of at 
least one of the visible final-state photons fulfil the following requirements: 

• > Pt° = 50 GeV (100 GeV) for at least one photon in case of W7 (Z7) produc- 
tion, 

• p l T > p l r£ = 20 GeV for all final-state charged leptons, 

• ARij > R c = 0.7, i,j = I, j, for the rapidity-azimuthal angle separation (|3.15|) 
between two charged leptons and between charged leptons and visible photons. 

• |?7| < rj c = 2.5 for at least one of the final-state photons and all final-state charged 
leptons. 

• p™ ms > j9™ lss,c = 50 GeV (100 GeV) for the missing transverse momentum in case of 
W7 (Z7) production for final states containing at least one neutrino, 

Furthermore, we introduce cuts on the invariant masses of the resonant gauge bosons. 
These cuts enhance the contributions of the gauge-boson pair-production process with 
respect to background diagrams and simultaneously reduce the error of the LPA. We 
impose the so-called reconstruction cut 

M z - 20 GeV < y/(p 3 +p 4 ) 2 < M z + 20 GeV (3.17) 
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on the phase-space for the process (|3.2|) . For the case of W7-pro duct ion (j3.3|) . the neutrino 
momentum is not an observable, such that one has to rely on the missing transverse 
momentum to restrict the squared energy of the lepton pair to the phase-space region 
around the resonance. To this end, a cut on the transverse mass 

M^<M w + 20GeV (3.18) 

is appropriate, where the transverse mass for a lepton-neutrino pair lv\ is defined as 

M% = ^/(K ss | + \pt,i\) 2 - (pf ss +p x ,i) 2 - (pf ss +Py,i) 2 - (3.19) 

For the case of the Z boson decaying into two neutrinos (J3.1|) none of the decay products 
is visible, and thus neither of the above reconstruction cuts can be imposed. 

For the bremsstrahlung process the scale Q is fixed using the possibly recombined 
momenta as follows. If both photons are visible, we use 

Q 2 = 1 {My + p 2 Ty + p 2 Tni + p 2 T J , V = Z, W, (3.20) 

to evaluate the PDFs. If only one photon is visible, the expression used for the scale 
Q depends on the considered final state. For the ZZ77 final state, we reconstruct the 
transverse momentum of the bremsstrahlung particle from momentum conservation and 
use the expression fl3.2()jl . When the final state contains one or more neutrinos, the 
bremsstrahlung momentum contributes to the missing momentum and we use f!3.13|) and 
dHHl) instead. 

3.2 Checks of the calculation 

To check the amplitudes generated by Pole numerically, we performed an explicit 
paper-and-pencil calculation of the Born amplitudes, the non-factorizable virtual correc- 
tions, and the real corrections. We have listed the corresponding expressions in App. [XJ 
In addition, we compared the Born amplitude and the bremsstrahlung amplitude for zero 
width numerically to the result of Madgraph [ [7TJ [72] for some phase-space points and 
found complete agreement. For a check of the amplitude for the factorizable virtual cor- 
rections, we generated the diagrams for the partonic processes corresponding to ()2.4)) by 
means of FeynArts. We then switched off the non-resonant and non-factorizable dia- 
grams by hand, and translated the corresponding amplitudes into Fortran code using 
FormCalc4. In the so-generated code, we replaced the resonant propagator by the cor- 
responding propagator including a finite width. The resulting amplitudes were checked 
by comparing to the output of Pole for some set of phase-space points, as well as by 
comparing the integrated cross sections, yielding in each case complete agreement. 

The numerical evaluation of the finite virtual corrections ()2.35|) was tested by numeri- 
cally varying the UV regulator, i.e. the mass parameter fi for the dimensional regulariza- 
tion, the infrared regulator, i.e. the infinitesimal photon mass m 7 , and the masses of the 
external fermions. The finite virtual corrections were found to be completely independent 
of these parameters, and we fixed the UV regulator to \i — 1 GeV, the infrared regulator 
to m 7 = 1 GeV, and the masses of the external fermions to the values listed in f!3.4|) . 
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Figure 1: Dependence of the sum of the singular virtual corrections, the complete real 
corrections, and the term arising from the MS definition of the PDFs for pp — > + X 
on the phase-space-slicing cut parameters. The horizontal lines indicate the error bars of 
our results obtained by using the subtraction method. 

To check the treatment of the real corrections, we calculated the dependence of the 
sum of ()2.16|) and (|2.17|) on the slicing cut parameters 5 S and S c and compared the cor- 
responding values to the sum of ()2.22|) and (J2.31J) obtained from the dipole subtraction 
method. In Figure Q we show the dependence on 5 S for fixed 5 C = 10~ 4 and the depen- 
dence on S c for fixed 5 S = 1CT 3 for the process pp — > v\l^ + X. The error bars reflect the 
uncertainty of the Monte Carlo integration. The cross section shows a plateau over the 
full considered range for 5 C , whereas the variation with 5 S leads to a plateau in the region 
10~ 5 < 5 S < 10~ 2 . In addition, the numerical values obtained in these ranges agree with 
those obtained by using the subtraction method within the integration error, which we 
indicate in Figure by the horizontal lines. The same analysis for the other considered 
processes yielded similar results. For the results to be presented in the next section, we 
accordingly fixed the slicing cuts to 5 B = 10~ 3 and 5 C = 10~ 4 . 

Finally, tree-level amplitudes and phase-space integration have been tested by com- 
paring the outcome of our code with the results obtained with the program used for the 
analyses in Ref. [|2H] • This latter code employs matrix elements generated by means of 
PHACT[F73|. a set of routines based on the helicity-amplitude formalism of Ref. [ I74j . and 
an independent integration method. For the sake of comparison, the original version of the 
program was extended to include the Z7-production processes (|3.1|) and (|3.2j) . The results 
for all three Born cross sections agree at the per-mille level. Furthermore, we numerically 
cross-checked the lowest-order results for the distributions to be presented in the next 
section for all three processes. Complete agreement has been reached, which confirms our 
implementation of phase-space integration, cut routines, and histogram generation. 
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4p } (fb) 


*& ),tot (fb) 


l/V^rg/ (%) 




212.26(7) 


-9.65(3) 


-4.5% 


0.5% 




38.99(9) 


-2.61(7) 


-6.7% 


1.1% 


Hi 


124.57(8) 


-2.36(5) 


-1.9% 


0.6% 



Table 1: Total lowest-order cross section (second column) as well as the electroweak 0(a) 
corrections in absolute size (third column) and in per cent of the lowest-order cross section 
(fourth column) for the three considered final states. The last column shows the statistical 
error for an integrated luminosity of L = 2 x 100/fb. 



4 Numerical results 

In this section, we illustrate the effect of the electroweak corrections on hadronic 
W7 and Z7 production at the LHC For the process (|3.1|) we sum over all three neutrino 
flavours, while for the two processes (|3.2|) and (|3.3|) we include the first two lepton families 
only. 

In the definition of the distributions for these processes, the possible presence of two 
photons in the final state, one coming from the real corrections, can lead to ambiguities. 
When observables of the final-state photon are involved, our prescription is to attribute 
the event weight to the bin corresponding to the photon with the largest energy. 

4.1 Total cross sections 

In Table [TJ we show the total hadronic cross section for the three processes (J3.1j) - (|3.3|) . 
The second column contains the lowest-order results, while the third and fourth entries 
display the value of the 0(a) corrections and their contribution relative to the lowest- 
order cross section, respectively. For the three considered processes, the electroweak 
corrections are negative and of the order of —2% to —7%. This is to be compared with 
the last column of Table ^ where we show an estimate of the statistical error based on 
an integrated luminosity L = 100 fb™ 1 for two experiments. As can be seen, the size of 
the electroweak corrections well exceeds the statistical uncertainty. They are of the same 
order of magnitude as the systematic error, which is expected to be dominated by the error 
on the PDFs and thus presently of the order of 5% for processes initiated by quarks. The 
systematic uncertainty could possibly be reduced by considering ratios of cross sections. 
On the other hand, the size of the O(a) corrections is strongly cut dependent and in 
fact increases with growing CM energy and vector-boson scattering angle. Larger values 
for the electroweak corrections are found by choosing, for instance, a more stringent cut 
on the transverse momentum of the photon. Such a constraint is useful for new-physics 
searches, although decreasing the statistics (see Ref. [|2H] and references therein). Also, 
distributions can be differently affected by the radiative corrections according to the 
specific observable at hand. In the following sections we illustrate this point for some 
sample variables. 
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Figure 2: Distributions in energy (left plot) and rapidity (right plot) of the final-state 
photon for the process pp — > v\Vi^ + X, I — e, /i, r. The inserted plots show the relative 
0(a) corrections in per cent normalized to the lowest order. 

4.2 Distributions for pp — > viVi^ + X 

We begin by considering distributions for pp — > v\Vfy + X. The error bars in the 
following plots result from the Monte Carlo integration errors. In the left plot of Figure |2l 
we show the energy distribution of the final-state photon. As previously mentioned, this 
distribution can be used to probe the existence of new particles which escape detection. 
These particles, once produced along with the photon, would lead to an excess of events 
in the high-energy domain. The kink at E 1 = 600 GeV is due to the cut \rj\ < on the 
rapidity of the final-state photon. This cut restricts the cosine of the photon production 
angle 9 1 in the laboratory frame to the interval 

exp- 2r?c -l exp 2r?c -l 

-2WCI1 <COS0 7 < W ~. r, 4.1 

exp zr ' +1 exp^ +1 

such that for our choice rj c = 2.5: 

^ = IsinfFl > 0.163071... « -. (4.2) 

Since we also applied a cut of p^ c = 100 GeV on the transverse momentum of the photon, 
this condition is always fulfilled for E 1 < 600 GeV and the rapidity cut is never applied for 
these photon energies. For E 1 > 600 GeV, on the other hand, the rapidity cut will start 
to exclude more and more events, which leads to a steeper decrease of the Born and NLO 
distributions. Assuming an LHC luminosity of L = 2 x 100/fb, the bin E 1 = 800 ±10 GeV 
collects a contribution of da/dE' 7 k5x 10~ 2 fb/GeV to the total cross section, and will 
thus contain about 

Hit 5x1 -2 fh 
N cv w L x — AEi = 200 ftT 1 x — 20 GeV = 200 (4.3) 
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events, where AE 1 = 20 GeV is the bin width. In this energy region, the electroweak 
corrections are of the order of —10%. The increase of the corrections with E" 7 can be at- 
tributed to large logarithms of Sudakov type. The electroweak corrections should therefore 
be included to match the experimental accuracy. This is the first calculation of the 0(a) 
corrections to the photon spectrum for viV\^j final states. Up to now, only QCD-corrected 
distributions have been computed. Since the electroweak corrections reduce cross sections 
and distributions, considering only lowest-order or QCD-corrected results could overes- 
timate the SM background when looking for new particles. An excess of events in the 
high-energy domain could then be misinterpreted as compatible with the SM predictions, 
and could therefore be missed. 

The right plot of Figure El shows the distribution in the photon rapidity. The elec- 
troweak corrections are here of the order of —3% to —5%, and are thus smaller as in the 
case discussed above and of the order of the present systematic uncertainty from PDFs. 
This is because all bins receive the dominant contributions from events with low CM 
energies, where the electroweak corrections are small. An estimate of the event rate as in 
f!4.3|) applied to the bin rp = ± 0.05 gives 900 events. 

4.3 Distributions for pp — > lly + X 

We turn now to the process (|3.2j) . mediated by Z7 production, and we show in Figure E] 
some sample distributions. The upper left plot displays the photon transverse momentum 
distribution. The kink at p\ = 250 GeV can be explained by the cut Ai? H - > R c on the 
rapidity-azimuthal angle separation of the two leptons resulting from the Z boson. The 
position of the kink is approximately given by 

pi = p% = = . 4.4 

v/2exp(i? c ) -1-1 

For p\ = p\ < pj? the R c cut is not effective, while for p\ = p\ > p^ the decay 
products of the Z boson are so strongly boosted that this cut eliminates many of the 
events. The effects of the electroweak corrections to this distribution have already been 
studied in Ref. [EU] for an on-shell final-state Z boson, and with the real corrections 
included in the soft and collinear limits only. In this case, the corrections were found to 
be large and negative, of the order of —20% and to increase in absolute size with rising 
p^. As can be seen from Figure El our calculation largely confirms these results. Since 
the effects of the anomalous couplings are most pronounced in the high-px region, the 
electroweak corrections could modify the experimental sensitivity to possible new physics. 
For the process considered here, the QCD corrections to the distribution in the photon 
transverse momentum have been calculated in Refs. [EHEEI^- ^ a PPly m g a jet veto, 
they are found to range between +60% and —20% [EJ. At NLO, the contributions arising 
from electroweak and QCD corrections are thus comparable in magnitude. Including the 
electroweak radiative effects in the experimental analysis is therefore mandatory in order 
not to overestimate the SM background. 

In the upper right plot of Figure El we show the distribution in the final-state invariant 
mass. A simple estimate as in ()4.3|) shows that the process should be observable up to 
= 750-1000 GeV. As can be seen from the inset plot, the total electroweak correc- 
tions amount to roughly —8%, independently of the value. Electroweak radiative- 
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Figure 3: Distributions in the transverse momentum of the final-state photon (upper 
left plot), in the 7Z invariant mass (upper right plot), in the photon rapidity (lower left 
plot), and in the photon-Z-boson rapidity difference (lower right plot) for pp — > + X, 
I — e, /i. The inserted plots show the relative G(a) corrections in per cent normalized to 
the lowest order. 
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correction effects are thus relevant for this important distribution. Earlier results for the 
electroweak corrections to the invariant-mass distribution have been published in Ref . [ ED] 
under the aforementioned approximations. Large negative corrections of the order of 
—15% were found, increasing in magnitude towards larger invariant masses. Comparing 
these results to the curve shown in the inset of the upper right plot in Figure El we rather 
find the relative size of the virtual and real corrections approaching a plateau with ris- 
ing invariant mass. We attribute this effect to angular-dependent logarithms, which are 
enhanced for events with photons collinear to the beams and compensate the Sudakov 
logarithms. Such events are eliminated by a cut on and thus do not influence the 

distribution, but contribute to the invariant-mass distribution also for large invariant 
masses. Thus, compared to the results published in Ref. [ED], the inclusion of the real 
corrections reduces the impact of the electroweak corrections for large invariant masses. 

The distribution in the photon rapidity is shown in the lower left plot of Figure 01 
The relative size of the electroweak corrections can be seen from the inserted plot to be of 
the order of —7%. Finally, the photon-Z-boson rapidity difference is shown in the lower 
right plot of Figure El While the corresponding corrections are at a similar level as for 
the two previous distributions in the central region, they become large and positive near 
|Arf z | « 4. For |Arf z | > 5, the cut | ?7 7 | < r] c = 2.5 eliminates all events with only one 
photon and a pair of leptons in the final state, and only bremsstrahlung events contribute. 
Thus, the smallness of the lowest order cross section causes the large relative corrections 
near and above |A?7 7Z | « 4. However, the event rate expected at the LHC in this region 
is too small for experimental observation. The remaining central part of the curve does 
not present novelties compared to the two previous distributions. 

4.4 Distributions for pp — > vil~f + X 

Turning next to the W7 production process, we show in Figure 0] some distributions 
for pp — > uil'-f + X, I — e, p,. Starting with the transverse-momentum distribution of the 
photon in the upper left plot, an estimate as in ()4.3|) shows that the event rate should 
be observable up to p^ = 300-400 GeV. An analysis of the effects of the electroweak 
corrections on this distribution summed over V[l and lv\ final states has already been 
published in Ref. [I2*3> taking into account the W-boson decay but including only the 
leading logarithmic virtual corrections to the production subprocess in the high-energy 
approximation. Large negative corrections were found, which can reach up to —12% 
at p\ = 400 GeV. As can be seen from the inset plot, we reproduce this result. At 
p^ = 400 GeV, the electroweak corrections are found to contribute about —13% of the 
Born result, which is well above the systematic and statistical errors. The QCD corrections 
to the transverse- momentum distribution of the photon have been analysed in Ref. [|2*T]. 
If applying a jet veto, large positive corrections of up to 60% for low transverse momenta 
have been found. They decrease to —5% to —10% of the Born cross section at p\ = 
400 GeV. For high transverse momenta, the electroweak corrections are thus of the same 
order of magnitude or even larger than the QCD corrections, and further reduce the lowest- 
order distribution. Again, a data analysis performed just including QCD corrections could 
overestimate the SM background in anomalous coupling measurements. 

In the upper right plot of Figure El we display our results for the distribution in the 
transverse mass of the final-state lepton-neutrino pair as defined in (J3.19|) . The curves 
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Figure 4: Distributions in the transverse momentum of the final-state photon (upper left 
plot), in the transverse mass of the lepton-neutrino pair (upper right plot), in the photon 
rapidity (lower left plot) and in the lepton-photon rapidity difference (lower right plot) 
for pp — > Z//Z7 + X, I = e, fi. The inserted plots show the relative electroweak corrections 
in per cent normalized to the lowest order. 
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terminate at the cut value M^ c = Mw + 20 GeV ~ 100.4 GeV. The W-boson resonance 
can be observed as a pronounced peak, an effect which could be used to determine the 
W-boson mass as a consistency check. The electroweak corrections can be seen from the 
inserted plot to contribute most at the peak, where they amount to —4% of the Born 
cross section. 

An interesting property of angular distributions for W7 pair production is the radiation 
zero, i.e. a kinematical configuration where all tree-level helicity amplitudes of the parton 
process are exactly zero. The radiation zero appears at 

cosd n 7T = q' ( 4 " 5 ) 

where 9 is the scattering angle of the photon in the partonic CM frame, and Q u and Qd 
denote the charge of the up quark and down quark, respectively. At the LHC, the radiation 
zero should be observable as a dip in the distribution in the rapidity difference between 
the photon and the charged lepton coming from the W boson [Ell^. It arises from gauge 
cancellations, and is characteristic for the SM. As a consequence, deviations from the SM 
gauge structure generally tend to fill in the dip, such that an analysis of the radiation 
zero provides an excellent opportunity to probe new physics. QCD corrections generally 
enhance the cross section at the radiation zero Depending on the process definition, 
the effect can be quite dramatic. If no jet veto is applied, the dip arising from the radiation 
zero is completely washed out. If selecting only events without final-state jets, a dip at 
the level of 20% survives. In the bottom right plot of Figure EJ we display our results 
for the distribution in the rapidity difference between the photon and the charged lepton, 
which exhibits the dip originating from the radiation zero at A7/ 7 « 0. The electroweak 
corrections hardly influence the dip and amount to about —2%. This is probably too small 
to have an observable impact on the data analysis. The smallness of the 0(a) correction 
has again its origin in the fact that a given bin receives the main contributions from the 
low-energy domain of phase-space, where the electroweak corrections are small. Selecting 
kinematical regions characterized by larger energies and scattering angles would increase 
the size of the electroweak corrections, at the same time enhancing the radiation zero dip, 
as shown in Ref. [12*3, but decrease the statistics. 

Analogous remarks apply to the electroweak corrections to the photon-rapidity dis- 
tribution shown in the bottom left of Figure |3J The central peak in this distribution is 
caused by events with small and disappears if the cut p^' c = 50 GeV is increased. 
The peaks in forward and backward direction are also present for larger p^' . They result 
from the peaking of the partonic matrix element in the direction of the incoming up quark 
which is additionally amplified by the boost from the CM to the laboratory system, which 
goes preferably along the direction of the valence quark. 

4.5 Comparison with virtual correction in high-energy approximation 

It is interesting to compare our results for the electroweak corrections to hadronic 
W7 production with those obtained in Ref. [EH]- I n this latter calculation, the virtual 
corrections were determined for the production subprocess only, and using the High- 
Energy Approximation (HEA) of Refs. [ [75J EUJ [77j . In this approximation, all terms that 
vanish with M^/y/I at high energies are neglected, and only the logarithmic contributions 
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Ref. [EH| 


this work 


^ c (GeV) 


4°p(fb) 

fir V / 


(T„ „ £*_*+._ (ft)l 
pp,nmtc v / 


rvirt 
pp. finite 


4°p(fb) 

t J t J \ / 


pp,nnitc v / 


rvirt 
pp, finite 


250 


6.01 


-0.292 


-4.86% 


6.02(3) 


-0.649(3) 


-10.8% 


450 


0.712 


-9.78 x lO- 2 


-13.74% 


0.711(4) 


-0.1372(7) 


-19.3% 


700 


9.30 x 10~ 2 


-2.06 x 10~ 2 


-22.12% 


9.31(6) x 10~ 2 


-2.54(2) x 10~ 2 


-27.3% 


1000 


1.25 x 10~ 2 


-3.72 x 10~ 3 


-29.6% 


1.253(9) x 10~ 2 


-4.32(3) x 10~ 3 


-34.5% 



Table 2: Lowest-order cross section and finite virtual corrections to the production sub- 
process for pp — > W 7 + X — > Ivfy + X, where lv% = lv\ + lv\ and I = e, \x. 

to the loop integrals of the form ln 2 (s/M^), ln(s/M^) hx(s/x) and hx(s/M^) are taken 
into account, where x = t,u denotes one of the Mandelstam variables for the production 
subprocess. The use of the HEA is justified by requiring a large transverse momentum 
of the final-state photon, which results in large values for the partonic CM energy \/§ 
compared to the W-boson mass, and thus to large values for the above logarithms. 

In order to tune our results to those of Ref. [EH]> we generated the amplitudes for the 
virtual corrections using Pole, but we switched off the diagrams involving the corrections 
to the decay subprocess. We evaluated the so-obtained amplitudes numerically by adapt- 
ing the input to the parameter and cut values used in Ref. [ EH] • For this comparison, we 
moreover summed up the process ()3.3|) and its charge-conjugate, i.e. we considered both 
W + 7 and W~7 production. The resulting numerical values for the Born cross section 
and the finite virtual corrections are shown in Table El for some values of the cut on the 
transverse momentum of the final-state photon. The numbers in parentheses represent the 
integration errors in the last digits. Note that the values shown in Table El do not exactly 
correspond to the results presented in a similar table in Ref. [EH]- This is because we 
found an error in the implementation of the azimuthal-pseudorapidity separation ()3.15|) 
when comparing our results to those of Ref. [EH]- The results for the Born cross section 
agree within errors at the per-mille level. Comparing our results for the finite virtual 
corrections to those obtained in the HEA in Ref. [EH], one would expect a sizeable dif- 
ference for low p^ c values, which decreases for higher values of this cut to the level of 
1-2%. However, the numbers in Table El point rather towards a constant difference of 
approximately 5-6%. 

In order to understand this effect, we used the fact that Pole employs the FF-based 
package LoopTools [[7S1IZH]; which numerically reduces all tensorial loop integrals onto 
a set of scalar basis integrals. We replaced each loop integral in this basis set by the 
corresponding expression in the HEA, and split the latter into different contributions. 
The total electroweak corrections obtained in this way are smaller than those in Table El 
by at most half a per cent. The sub contributions obtained with this procedure are shown 
in Table El 

The electroweak logarithmic contributions in the second and third column of Table El 
result from the contributions to the basic loop integrals that involve only logarithms of the 
form ln 2 (s/M^), ln(s/M^) ln(s/x), x = i,u and ln(s/M(y). This is the approximation 
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p^' c (GeV) 


a EWLog(fk) 


°EWLog 


AagUfb) 


°QED 


O-rcm(fb) 


r(l) 

L'rem 


250 


-0.364(1) 


-6.05% 


8.6(1) x 10~ 2 


1.43% 


-0.360(2) 


-5.99% 


450 


-0.1026(4) 


-14.4% 


1.13(2) x 10~ 2 


1.59% 


-4.29(2) x 10~ 2 


-6.04% 


700 


-2.11(1) x 10~ 2 


-22.6% 


1.40(3) x 10~ 3 


1.50% 


-5.31(3) x 10~ 3 


-5.71% 


1000 


-3.77(3) x 10~ 3 


-30.1% 


1.74(5) x 10~ 4 


1.39% 


-6.63(5) x 10~ 4 


-5.29% 



Table 3: Finite virtual corrections to the production subprocess for pp — > W ± 7 + X — > 
Ivf-j + X in HEA split into electroweak logarithmic contributions 0EWLog' the QED rem- 
nants A<Tqe D , as well as the remaining contributions arlL- 

used in Ref. [EHj. As can be seen by comparing these numbers to the results in the third 
and fourth column in Table El the leading logarithmic contributions differ by about 1.2% 
for = 250 GeV and this difference decreases to the per-mille level for higher cut values. 
Thus we can reproduce the results of Ref. [|2H] to a satisfying accuracy by keeping only 
the electroweak logarithms. 

The QED contributions in the fourth and fifth column of Table El are obtained by 
keeping only terms involving logarithms in the regulators for the masses of the photon 
and the external light fermions, i.e. logarithms of the form m(m^/M{^) or ln(m^/M^), 
in the loop integrals and subtracting the finite virtual corrections as given in ()2.39|) . We 
here observe a roughly energy-independent contribution of about 1.5%. This contribution 
results from the fact that the subtraction term ()2.39|) involves logarithms of the form 
ln(m^/sjj) and \n.{mf/ Sij) while in the HEA they appear as ln(m^/M{y) and hi{mf/M^j). 

Finally, the contributions listed in the sixth and seventh column of Table El result from 
all remaining terms in the HEA, i.e. logarithms depending only on ratios of kinematical 
variables and constant terms in the high-energy expressions of the basic loop integrals, as 
well as extra constant terms arising from the reduction of tensor loop integral to scalar 
loop integrals. These remaining terms yield a contribution of 5-6% and thus make up the 
largest part of the difference between our results and those of Ref. [ • 

In order to trace the origin of this large effect, we further split the remaining con- 
tributions to the basic loop integrals as shown in Table EJ The numbers in the fourth 
and fifth column comprise the contributions linear in the purely angular-dependent log- 
arithms, i.e. logarithms of the form ln(s/x), x = t,u, whereas the numbers in the sixth 
and seventh column correspond to the terms containing these logarithms squared. As 
can be seen from Table EJ the contributions linear and quadratic in the purely angular- 
dependent logarithms vary only little with increasing energy, but contribute with about 
—2% and —(3-4)%, respectively, i.e. the make up the bulk of the remaining contribu- 
tions. The large size of the contributions of the purely angular-dependent logarithms can 
be explained by the enhancement of the cross section of the considered processes in the 
forward and backward directions. Since in these regions s/t or s/u are small the purely 
angular-dependent logarithms become large. The contributions of these logarithms are 
the source of comparably large difference between the corrections in the HEA of Ref. [ EH] 
and the results of this paper. On the other hand, as shown in the second and third column 
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p^(GeV) 


^rem.NoLoglfk) 


r(i) 

rem,NoLog 


^rem.Loglfk) 


6 m 

rem, Log 


^.Log^fb) 


U T 1 

rem Loff 


250 


-1.55(8) X 1CT 2 


-0.26% 


-0.107(1) 


-1.78% 


-0.238(2) 


-3.95% 


450 


-4.13(8) x 10~ 3 


-0.58% 


1.139(2) x 10~ 2 


-1.97% 


-2.48(2) x 10~ 2 


-3.50% 


700 


-6.94(9) x 10~ 4 


-0.74% 


-1.78(2) x 10" 3 


-1.91% 


-2.83(2) x 10~ 3 


-3.05% 


1000 


-1.09(1) x 1CT 4 


-0.87% 


-2.25(3) x 10~ 4 


-1.79% 


-3.29(3) x 10~ 4 


-2.62% 



Table 4: The non-leading contributions to the virtual corrections to the production sub- 
process for pp — > W ± 7+X — > lvi^+X in HEA split up into terms cr^ NoLog containing no 

angular-dependent logarithms, into terms cr^ m Log depending on single angular-dependent 
logarithms as well as into terms ^emLog 2 depending on angular-dependent logarithms 
squared. 

of Table 0] the non- logarithmic remaining terms contribute to the cross section only below 
the per-cent level. The increase of this contribution with rising transverse momentum cut 
of the photon can be attributed to the energy dependence of those parts of the matrix 
elements that have not been evaluated in the HEA. 

We thus find that, for processes dominated by contributions coming from phase-space 
regions with small kinematical variables, the accuracy of the HEA which only takes into 
account enhanced logarithms of s/M 2 may only be at the level of several per cent, unless 
the enhanced angular-dependent logarithms are included as well. 

5 Summary and Conclusions 

We have calculated the electroweak corrections to Z7 and W7 production at the LHC, 
with Z and W decaying into leptons. For the Born cross sections and the real corrections 
we use complete matrix elements. The virtual corrections are included in the leading-pole 
approximation. For a treatment of the soft and collinear divergences, we used the dipole 
subtraction method as well as the phase-space-slicing technique. 

We have implemented our strategy in a Mathematica package called Pole, which works 
as an extension of the programs FeynArts3 and FormCalc3 . 1 and is designed to calculate 
electroweak corrections in leading-pole approximation for hadronic or partonic processes. 
We have applied the program to the pp — > v\Vfi , ll'y , ujj processes at the LHC, and 
tested the results extensively. As part of our checks we have performed a comparison 
with the virtual corrections evaluated in the high-energy approximation, which reveals 
that angular-dependent logarithms can give effects of several per cent at high energies. 

For the above-mentioned processes and typical LHC cuts, we find electroweak correc- 
tions of the order of —5% for total cross sections and angular distributions as, for instance, 
the radiation zero in W7 production. For transverse-momentum and energy distributions, 
the 0(a) corrections contribute up to —20% of the lowest-order cross section. Their in- 
fluence is thus well above the systematic and statistical errors, and may be of the same 
size as the QCD corrections or even larger, depending on the cut settings. Consequently, 
electroweak radiative effects have to be included in the experimental analysis when look- 
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Figure 5: Tree- level diagrams for q\ + q 2 — > f% + /4 + 7 

ing for the existence of anomalous vector-boson couplings or for the production of new 
particles via hadronic production of Z7 and W7. Thus, the Monte Carlo generators for 
LHC should include besides QCD corrections also the electroweak corrections. 
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A Analytical results for the amplitudes 

In this appendix, we summarize analytical results for those parts of the amplitudes that 
are compact enough to be explicitly displayed. We present the analytical expressions for 
the Born amplitude in Section IA.11 the amplitude corresponding to the non-factorizable 
corrections in Section IA.21 and the amplitudes for the QED bremsstrahlung processes 
in Section TA. 31 The calculation of the Born and bremsstrahlung amplitudes are carried 
out in the framework of the Weyl-van der Waerden formalism [ IHUJ, as presented in 
Ref. [E2|. 

A.l The Lowest-order amplitudes 

The helicity amplitudes of all partonic processes (|2.4|) can be obtained from the generic 
set of Feynman diagrams shown in Figure 03 The amplitude corresponding to these dia- 
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Table 5: Coupling constants of quarks and leptons to the gauge bosons V = 7, Z, W ± 
grams can be written as 

E Cv\ 2qi c;% f Br 3X5 (PuP2,P3,P4,Ps, Q qi ,Q q2 , QfrQu), (A.i) 

V 

where the quantities Qi,i = 1, ...4 denote the charges of the external fermions (the 
charges of the external antifermions are given by —Qi.) and ci,C2 are the colour indices 
for the initial-state quarks. The sum over V runs over all SM vector bosons. The non- 
vanishing gauge couplings Cir? of a fermion fa and an anti-fermion f a to the vector boson 

* Jajb 

V are listed in Table El using the following conventions. The particle indices I (q) denote 
leptons (quarks). The charges and the third component of the isospin of the fermion / 
are denoted by Qf and Ij, respectively. The quantities s w and c w are the sine and the 
cosine of the weak mixing angle, whereas a is the fine-structure constant. Finally, 5ij and 
Vij denote the Kronecker delta and the quark mixing matrix. Note, that the elementary 
charge e has been extracted from the couplings. 

As discussed in Section 1231 the amplitude can be split into a resonant and a non- 
resonant part. This splitting can be done on the level of the generic functions 

go-i<T3A 5 grcs,o-i<T3A 5 _|_ gnres,<Tio-3A 5 ^ 2^ 

The resonant part comprises all contributions from diagrams containing the resonant 
propagator of the decaying particle V, and thus corresponds to the amplitude originating 
from the diagrams in the first row of Figure El Applying the Weyl-van der Waerden 
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formalism as presented in Ref . [ |82| to these diagrams, the generic function for all helicities 
positive is given by 

{PiPa) 2 



By"' (Pl,P2,P3,P4,P5,<9l,Q2,<23,<34) = -Pv(P3 + Pi) 



(pm) 



x {%<^ + %<^-(Q 2 -Q 1 )i'v( Pl + P2 )(p 3 p 5 >*( P2Pl )*V (A.3) 

The non-resonant part of the Born amplitude, on the other hand, receives contributions 
from the remaining diagrams in the second row of Figure El Expressed in Weyl-spinor 
products, the generic function for all helicities positive reads 

D nrcs,+++/ r\ r\ r\ \ D ( , \Qs(PiP4} 2 (P2Pi}* { A A\ 

By \Pi,P2,P3,Pa,P5,Qi,Q2,Q3,Qa) = PviPi + P2)—, r~, r — • (A.4) 

{P3P5)(P4P5) 

Here My is the mass and IV the decay width of the vector boson V, and the corresponding 
propagators are denoted as 

= *»-«?,-IM,IV (A ' 5) 



The Weyl-spinor products are defined by 



(pq) = 



t>„ dp ■ dq _y, Oq . Op 

p cos — sm — — e Vq cos — sm — 
2 2 2 2 



(A.6) 



where 9 P , 9 q and 4> p , <p q are the polar and azimuthal angles, respectively, of the corre- 
sponding light-like 4-momenta 

P 11 — Po(l) sin 9 P cos <j> p , sin 9 P sin </> p , cos 9 P ) , 

9 M = qo(l, sin 9 q cos (fi g , sin 9 q sin <fi g , cos 9 q ). (A. 7) 

The amplitudes for neutral current reactions can be extracted from (jA.3|) and (jA.4|) by 
simply setting (Q2 ~Qi) — 0. Inserting ()A.3|) and (|A.4|) in ()A.1|) and adding the resulting 
expressions yields the full Born amplitude. While this is gauge invariant, the resonant 
and non-resonant parts as defined above are gauge-dependent, unless applying the LPA. 
To this end, one has to set the finite width in the s-channel propagator Py{pi +P2) in the 
resonant amplitude (jA.3|) to zero and to insert the resulting expression in (jA.lj) . where the 
spinor and scalar products have to be evaluated using the on-shell projected momenta, 
except for in the resonant propagator Py{p% +P4). Note that we set the gauge spinor for 
the external photon to g = to obtain the above results, which is why the contribution 
of the diagram with the photon coupling to the external fermion / 4 vanishes. 

Dropping the superscript for resonant and non-resonant contributions in what follows, 
the other helicity combinations can be obtained as 

By ++ (pi,P2,P3,P4,P5,Ql,Q2,<33,<24) 

= B v ++ (p 2 ,Pl,P3,P4iP5, -Q2, -Ql,Q3,Ql), 
By" + (pi,P2,P3,P4,P5,Ql,Q2,Q3,<24) 

= B y ++ (pi,p 2 ,P4,P3,P5,<2l,<22, -Q4, -Q3), 
By~ + (pi,P2,P3,P4,P5,Ql,Q2,Q3,<?4) 

= B v ++ (p 2l p u p 4 ,p 3 ,p 5 , -Q 2 , -Q u -Q 4 , -Q 3 ), (A.8) 
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as well as 



By 1<T3 (P1,P2,P3,P4,P5, Ql, Q2, Q3, Qi) 



[B v ai CT3+ (pi, p 2 , P3, Pa, Pb, Qi, Q2, Qs, Q* 



(A.9) 



A. 2 The non-factorizable corrections 

A general discussion for the photonic non-factorizable corrections for a process involv- 
ing massless external particles and an arbitrary number of resonances can be found in 
Ref. [Ell- Applying the general results displayed there to our processes, the correction 
factor results in 



5^(n s \(p 3 +p 4 ) 2 -M v ) 



E T * T jQiQj- Re { AV (pT h ; (pa + PA),pf h )} , v = w,z (a.io) 

t=l,2 i=3,4 
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A v (p? h ;(p 3 + Pi ),p° sh ) 



2 In 



Ml -iM v T v -(p 3 +p 4 y 



ln|-|-l 



1- J ) 



The kinematical variables i and t are defined by 



* = (pT h - K sh ) 2 ' * = (pT h - p? h - p°tf - M? 



v ■ 



(A.12) 



and Li 2 denotes the dilogarithm (|2.21j) . 



A. 3 The real QED corrections 

The determination of the helicity amplitudes for the real QED corrections can be 
performed in the same way as described in Section lA.ll for the case of the Born amplitudes. 
The calculation is much more involved, though, since the set of generic diagrams comprises 
here 31 instead of five diagrams. These diagrams can be obtained from the ones shown in 
Figure in the following way. A total number of thirty diagrams is obtained by attaching 
the additional real photon to each fermion line and to each massive vector-boson line 
in the Born graphs shown in Figure The one remaining diagram results from both 
photons coupling to the intermediate vector boson via a VV'yy coupling. Apart from the 
increased number of diagrams, one has to calculate the generic amplitude for two different 
polarization combinations, since a final state with both photons having the same helicity 
is not related by any discrete symmetry to a final state with two photons of opposite 
helicity. 

Adopting the conventions of Section IA.1I the generic amplitude of the QED brems- 
strahlung process reads 



MZ^dTcLc, (Pi. Pa, P», P4, P5, Ps) 



-4e 5 S S 



E 



^ 29l ^> 3/4 By i ; 3A5A6 (pi,P2,P3,P4,P5,P6, Q qi ,Q q2 ,Q h ,Q h ] 



(A.13) 
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where the sum over V is again taken over all SM vector bosons. The couplings of the 
gauge boson V to the fermion f a and the anti-fermion f b is again denoted by C'' f , and 

* lb J a 

can be read off Table |5j The generic function for all helicities positive is obtained as 



B v^ ++ (Pi ' V2 , P3 , Pa, P5 , Pe , Qi , Q2 , Q3 , Qa) 



(P4Pl) S 



x { Pv(Ps + Pa) 



{PAPb){PAP&) 

-QlipmY : {paPi) , QiQ2{pa\p2 + psIps) 



{pm)(piPs) (P2P$){pm) 

QiQ2{Pa\P2+P0\P3) , Ql{PlP-iY (P2Pa) 



(PiP5)(P2Pg) 



(PsPs) 



(P2P5)(P2Pe) 



P V (P3+Pa)(Q2-Qi)(P5P3Y 



Qi(Pa\pi +Pe|P2> Q2(P4|P2 + Pe|pi) 



(PlPe) 



(P2Pe) 



^V(Pl + P2XQ2 - Ql)(PlP2>*(P4|P2 +Pl|P6> 



+ iV(P3 + P4 +Pe) 



<2a 



(PsPe) 



+ Py(P3 + p 4 )(Q 2 - Ql){P6P3Y 



Qi(Pa\Pi+P5\P2)_ Q2(pa\P2 +Ps\pi) 



(pipb) 



(P2P5) 



Pv(Pl+ P2)(Q2 ~ Ql){PlP2Y{Pi\P2+Pl\Pb) 



Pv{Pl+P2)- 



■(PiP2Y(P3Pa) 



(P3Pb) (psPe) 

where in the above expressions, a product like {pi\pa + P&\Pb) is defined by 

(a\b + c\d) = (ab) (db)* + (ac) (dc)*. 

The corresponding function for helicity A5 negative reads 

1 

B £f + (Pl, P2, P3, P4, P5, P6, Ql, Q2, ^3, Q4) 



X < P V (P3 +Pa) 



(Pm)(P3P5) 

Qi {P3P2Y 



(A.14) 



(A.15) 



+ 



QI 



ipmYipm) (pi +P5 + pg) 2 

f" (PlP4 ) (P1P5 ) * (PlP5 ) (P2P4 ) (P3^2 ) 

{PaPi} 



(Pa \pi + P5 IPs) (Pa \pi + P5 be) (piPe) 



(P2PbY{P2Pf>) (P2 +P5+ Ps)' 



(P5\P2+P(>\P3}(Pa\P2+P6\P3}(P2P5Y 



33 



+ (PlP3 ) * (P2P&Y (P2P6 ) (P2P3 ) * (PiPl ) 



+ ~ -^tt^ r(pm) 2 {{P2P?>Yf - ~, r{Pi\Pi +Psbs) 2 



(P2P5}*{P1P6) 
+ Pv(P3 + Pi + P5){Pl\P4 + P5IP3) 



{PiPsY 



{PlPb)*{P2P&) 



Pv(P3+Pi)(Q2 ~ Ql){PbPi) 



X 



Q 



(PiPa) 



1 -t(P4Pi)(P3P2Y 



Qi 



(P2P6} 



(Pi\p2 + Pa\p3) 



+ p v (pi +P2YQ2 - Qi) (pmY (pePsY (pm) + (p*\p2 +pi\p%)(pmY 



+ Pv(P3 +Pi+ Pe)(P4|p3 + Pe|p2> 



(PsPe) 



Pv(P3+P4)(Q2-Ql)(P6P3Y 



X 



/ ® 1 x * {p4p2+pg\pz) + 1 ® 2 ^ (pm)(P3P2Y 
(P1P5) \P2P5Y 



Pv{Pl+P2){Q2 - Ql)y(P5Pl)(P5P4)(P5P3Y + (P5\P2 + Pl\Ps) (PiPl) 



+ P V {PI+P2) 



Pv(P3+Pi)(Q2 ~ Qi) 2 {PiP4,)(P2P3)*(P4P5)(P6P3Y 



Ql (PdPi )(P3P6 YiPiPl / I . I v 

+ I T / — 7 \P5 \P3 + P6 \P2) 

We) (P3+P5+P6) 2 

Q3Q4 



(Pi \pi + P5 \Ps) (Pi \P3 + P6 |P2> 

(pi\Pi + PslPe) 



{P3Ps){PaPs) 

Ql (P5P3Y 'ipmY : (pm) 



{PiPbY (Pi+Ps+Pe) 2 
The other helicity combinations can be calculated by 

r + 



(A.16) 



B y l" A5+ (pi, p 2 , P3, Pi, P5, Pg, Qi, Q2, Q3, Qi) 



By+ A5+ (p 2 ,Pl,P3,P4,P5,P6, -Q2, -Ql,Q3,Q4), 



B V 7 A5+ (Pl, P2, P3, P4, P5, P&, Qli Q2i Q3, Qi) 



By^ A5+ (pi,P2,P4,P3,P5,P6,Ql,Q2, ~Qi, -Q3), 



By 7 5 (pi,P2,Pa,Pi,P5,Pe,Qi,Q2,Qa,Qi 

V7 



By+ A5+ (p 2 ,Pl,P4,P3,P5,P6, -Q2, -Ql, -Qi, ~Q- 



(A.17) 



and 



By 7 ff3 A5 (Pl,P2,P3,Pi,Pb,P6,Ql,Q2,Q3,Qi) 
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[B^ 1 ff3 A5 + (pi,p 2 ,P3,P4,P5,P6,Ql,Q2,Q3,<54 



B Explicit form of the (Lorentz-invariant) on-shell projection 



. (A.18) 

P^Pv 



In order to be able to evaluate the virtual corrections to the processes ()2.3|) in LPA, 
an on-shell projection has to be specified which maps a general set of momenta $o on a 
set of momenta $Q Sh such that (p^ sh +pf h ) 2 = M v , where My is the mass of the decaying 
vector boson V = W, Z. Such a projection is given by 

p3 = Vm, P5 = V2P5, P4=Pl+P2-P3-Pb, (B.l) 

where scaling variables x and z are obtained from the mass-shell conditions for the reso- 
nance and for the momentum p 4 as 

S-My My 

y 2 = —, \ — , yi = —, : — . (B.2J 

2(pi + P2)Pb 2{Pi +P2- zp 5 )p 3 

This on-shell projection only involves Lorentz-invariant quantities. 
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